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Abstract. The decline of the high mass X-ray transient V 0332+53 during the Dec. 2004 to Feb. 2005 outburst is analysed 
from the data recorded by INTEGRAL. The flux is shown to decrease exponentially until 2005 Feb. 10, with a decay time scale 
of ~30 days above 20keV and ~20 days at lower energies, and to decrease linearly thereafter. 

The energy spectrum is well modelled throughout the decay by a power law with a folding energy of ~7.5 keV, and with two 
cyclotron absorption features. The folding energy does not vary significantly over the decay, but the spectrum becomes harder 
with time. Most importantly, we show that the parameters describing the fundamental cyclotron line around 27 keV do vary 
with time: its energy and depth increase (by about 17% for the energy in ~6 weeks), while its width decreases. These changes 
of the cyclotron line parameters are interpreted as resulting from a change in the extent of the cyclotron scattering region. 

Two quasi-periodic oscillations are also observed at various times during the observations, one at 0.05 Hz and another one near 
the pulsation frequency around 0.23 Hz. 
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1. Introduction 

The recurrent X-ray transient V 0332+53 experienced an X- 
ray outburst from December 2004 to February 2005. The out¬ 
burst was predicted one year earlier from an optical brightening 
of the massive 08-9Ve companion star, BQ Cam (Goranskij 
& Barsukova 1 and was first detected by the Rossi 
X-ray Timing Explorer ( RXTE ) in November 2005 (Swank 
et al. 2504;). The INTernational Gamma-Ray Astrophysics 
Laboratory ( INTEGRAL ) started to observe the source on 
2005 Jan. 6 (Kreykenbohm et al. [20051 hereafter Paper I) and 
monitored its X-ray activity until the end of February. 

The observation of this fourth detected outburst of 
V 0332+53 in thirty years revealed the presence of three cy¬ 
clotron lines (Coburn et al. 120051 Paper I, Pottschmidt et al. 
I2005li . The fundamental has an energy of 28keV, and the first 
and second harmonics are at 50 keV and 71 keV, respectively. 
The magnetic field strength is 2.7x 10 12 G (Paper I, Pottschmidt 
et a I. f2005k The spectrum is otherwise described by a power 
law with a high-energy cut-off, typical of the spectra of X-ray 
pulsars. A new Quasi Periodic Oscillation (QPO) around the 


0.23 Hz spin frequency of the neutron star has also been re¬ 
ported in the RXTE data (Qu et al. :2005j), in addition to the 
0.05 Hz QPO discovered by Takeshima et al. ( 119941 in the 
Ginga data of the 1989 outburst. 

In Paper I, we presented the INTEGRAL observation of 
V 0332+53 around the outburst peak. This paper presents the 
evolution of V 0332+53 during its decline phase. The obser¬ 
vations and data analysis method are described in Sects Q and 
[3 respectively. We then present in Sect.|4]the evolution during 
the decline of the fluxes in several energy bands, of the energy 
spectrum, and of the variability spectrum of V 0332+53. The 
results are then discussed in Sect.|^ 

2. Observations 

The INTEGRAL satellite carries three high energy instru¬ 
ments, the imager IBIS (operating in the range 15 keV to 10 
MeV), the spectrometer SPI (20 keV to 8 MeV) and two iden¬ 
tical X-ray instruments JEM-X1 and JEM-X2 (3 to 35 keV) op¬ 
erational in turn (see Winkler et al, 12003l for more information 
on the INTEGRAL mission and instruments). IBIS is further 
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Table 1 . Summary of observations. From first to last column: (1) Revolution and pointing numbers; (2) Earth date of the start 
of the observation in UTC / IJD (INTEGRAL Julian Day = MJD -51544) / MJD (Modified Julian Day); (3) Earth date of the 
end of the observation in MJD; (4) Individual pointing duration; (5) Total IBIS observation time; (6) Sum of IBIS Good Time 
Intervals; (7) Total JEM-X1 observation time; (8) Sum of JEM-X1 Good Time Intervals. 


Revolution / 
pointings 

Start date 

(UTC / IJD / MJD) 

End date 
(MJD) 

pointing 

duration 

(sec) 

IBIS 

Exposure 

(sec) 

X GTIs 
(sec) 

JEM-X 1 

Exposure X GTIs 
(sec) (sec) 

272 / 75 

2005-01 -06T08:48 / 1832.36724 / 53376.36724 

53376.76044 

42184 

42184 

30308 

36896 

36446 

273/66- 81 

2005-01-08T22:18 / 1834.92997 / 53378.92997 

53379.65749 

3500 

60953 

59586 

60953 

56693 

274 / 04 - 09 

2005-01 -10T03:20 / 1836.13991 / 53380.13991 

53380.39113 

3500 

21119 

20901 

21119 

19801 

278 / 48 - 74 

2005-01-23T16:39 / 1849.69471 / 53393.69471 

53394.66610 

2858 

80901 

58411 

75190 

73078 

284 / 04 - 45 

2005-02-09T01:15 / 1866.05280 / 53410.05280 

53410.86456 

3572 

149086 

147739 

149086 

147469 

285 / 02 - 08 

2005-02-12T01:02 / 1869.04402 / 53413.04402 

53413.23174 

2200 

15517 

15349 

15517 

14394 

286/ 18 - 24 

2005-02-15T10:41 / 1872.44588 / 53416.44588 

53416.63218 

2200 

15403 

15093 

15403 

15185 

287 / 02 - 08 

2005-02-18T00:37 / 1875.02612 / 53419.02612 

53419.23077 

2200 

16988 

16628 

16988 

15809 

288 / 02 - 08 

2005-02-21T00:25 / 1878.01820 / 53422.01820 

53422.25899 

2200 

20101 

19226 

20101 

18854 


composed of two modules, ISGRI operational down to 15 keV 
and PICsIT above 200 ke V. 

INTEGRAL observed V 0332+53 from 2005 Jan. 6 to Feb. 
21. Three main observation sets were scheduled, in revolutions 
272-274 (~120ksec), 278 (~80ksec) and 284 (~150ksec), fol¬ 
lowed by four ~15ksec observations from revolution 285 to 
288 to monitor the tail of the outburst decay. The observations 
are summarized in Table 0 Except for the observation in rev¬ 
olution 272, which was performed in a staring mode, all ob¬ 
servations were made following a hexagonal dithering pattern. 
The hexagonal dithering pattern consists of one pointing (sta¬ 
ble spacecraft attitude) with V 0332+53 positioned at the center 
of the held of view of the instruments, and six pointings with 
V 0332+53 off-axis in steps of 2 degrees in a hexagonal pat¬ 
tern. The integration time for each pointing may vary from one 
observation to the next and is given in Tabled 

The V 0332+53 observations were affected by solar activ¬ 
ity which impacted on the quality of the data. In particular, a 
strong solar hare lasting from Jan. 20 to 23 lead to a split of 
the IBIS data of revolution 278 into 7760 small Good Time 
Intervals (GTIs) of only few seconds each. The JEM-X GTIs 
were greater than 97% of the observation time for that revolu¬ 
tion, and in general greater than 90% for all revolutions (see 
Table Q. 

3. Data analysis 

We analyse all the public data collected by the ISGRI and the 
JEM-X 1 instruments on board of INTEGRAL. PICsIT is not 
useful for our analysis because the spectrum of V0332+53 falls 
off above lOOkeV, and SPI data are not used either since its 
sensitivity in the 20 to 100 keV range is less favorable than 
ISGRI (Paper I). Pointings close to perigee passage times are 
further discarded from the analysis as they are contaminated by 
the Earth radiation belt. 

The ISGRI and JEM-X data are processed with version 5.0 
of the Offline Scientific Analysis (OSA) package delivered by 
the INTEGRAL Science Data Center (http://isdc.unige.chi us¬ 
ing default settings. Images and lightcurves are produced in 


the 3-5 keV, 5-10 keV and 10-15 keV energy bands for JEM- 
X, and in the 20-30 keV, 30-40 keV and 40-60 keV bands for 
ISGRI. 

Contamination from other sources is not a problem as 
V 0332+53 is the only bright X-ray source in the ISGRI field 
of view. Systematic errors are however known to affect the 
data analysis. In order to estimate the amplitude of those sys¬ 
tematic errors, and to calibrate the results, an analysis of the 
Crab was performed with OSA parameters identical to those 
used for the analysis of V 0332+53. Crab public data are taken 
from revolution 239, pointings 3 to 27 (Modified Julian Day 
MJD=53275.5289 to 53276.1326), during which the Crab is lo¬ 
cated within 3 degrees from the center of the held of view. The 
resulting mean Crab fluxes extracted from the images in the 3- 
5, 5-10, 10-15 keV bands of JEM-X are 43.5 cps (14% error), 
43.2 cps (10%) and 16.4 cps (12%), respectively. The quoted 
errors represent the maximum deviations measured from the 
mean fluxes over the 25 pointings. The Crab fluxes in the 20- 
30, 30-40 and 40-60keV energy bands of ISGRI are 78.8 cps 
(13%), 39.4 cps (7%) and 45.9 cps (4%), respectively. 


4. Results 

4.1. Fluxes 

The fluxes of V 0332+53 in the ISGRI and JEM-X energy 
bands are plotted in Fig. |]]as a function of time. They are nor¬ 
malised to the Crab fluxes in the respective bands. 

The decline phase can be described by an exponential de¬ 
cay up to MJD^53412, followed by a linear decrease. The ex¬ 
ponential decay is less rapid at the ISGRI energies than at the 
JEM-X energies. From a simple estimation, the exponential de¬ 
cay time is about 30 days above 20 keV (solid curves in Fig. m 
and about 20 days below 15 keV (solid curves in Fig. m- The 
exponential decay is also confirmed by the 2-12 keV lightcurve 
recorded by the All Sky Monitor (ASM) on board of RXTE, 
plotted in Fig. Ub for comparison with the JEM-X lightcurves 
(see Levine et al. l$96t for more information on the RXTE satel¬ 
lite). 
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Fig. 1. a. Pointing averaged ISGRI flux evolution of V 0332+53 
as a function of time in the different energy bands as labelled in 
the figure. The fluxes have been normalised to Crab units. The 
error bars displayed in the figure represent only the statistical 
errors; neither do they take into account the uncertainties in the 
Crab fluxes in the respective energy bands. The solid curves are 
exponential (from MJD=53376 to 53412 with r=30 days) and 
linear (from MJD=53412 to 53422) decays superimposed on 
the data. b. Same as a, but for JEM-X fluxes and with an expo¬ 
nential decay of r=20 days for the solid curves. The 2-12 keV 
RXTE/ASM lightcurve is also shown, with a normalisation fac¬ 
tor of 77 cps for 1 Crab. 


The evolution of the ISGRI and JEM-X hardness ratios are 
shown in Fig. |3 Above 20keV the spectrum of V 0332+53 
softens with time, with the [30—40] keV / [20-30] keV ratio de¬ 
creasing by about 36% from an average value of ~0.55 in the 
first observation to -0.35 in the last observations. On the other 
hand, below 15keV, the spectrum becomes harder, with the 
[10-15] keV /[5-10]keV ratio increasing by about 28% from 
an average value of ~ 1.8 to -2.3 during the decline. As will be 
seen in Sect. 14.21 the opposite behaviour of the hardness ratios 
above 20keV and below 15 keV are understood from the evo¬ 
lution of the spectral shape of V 0332+53 during the decline. 

The ISGRI/JEM-X hardness ratio ([20-30] keV/[5- 
15] keV) also evolves with time, increasing from an average of 
-0.7 to -1.5 during the decline. This results from the increase 
of the spectral power law index presented in the next section. 
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Fig. 2. Hardness ratios computed from the pointing-averaged 
(and Crab normalized) fluxes displayed in Fig. [I] Upper trian¬ 
gles: ISGRI [30-40 keV]/[20-30keV]; lower triangles: JEM-X 
[10-15keV]/[5-10keV]; squares: ISGRI [20-30keV] / JEM- 
X [5-10keV]. 


4.2. Energy spectra 

The JEM-X and ISGRI spectra of V 0332+53, normalised to 
the Crab spectrum, are shown in Fig. 0 at four different times 
during the decline of the source. The maximum flux relative to 
the Crab is recorded between 10 and 20keV, at an energy E max 
which increases with time. Fitting gaussian curves to the JEM- 
X spectra in a lOkeV width range around the maximum leads 
to Emax = 14.6, 15.4, 16.8 and 18.0keV for revolutions 273, 
278, 284 and 287, respectively. This shift of E max explains the 
hardening (softening) at low (high) energies noticed from the 
JEM-X (ISGRI) hardness ratios presented in Sect. 14. II 

The spectra of accreting X-ray pulsars include contribu¬ 
tions from at least an accretion column and from the hotspots 
at the surface of the neutron star. In the absence of a consis¬ 
tent model, various empirical models are commonly used in 
the literature to fit such spectra (see for example Kreykenbohm 
et al. 119991 for a discussion). In this study, we use the simple 
“cutoffpi” model of the form 

F(E) = A • E~ r ■ e ~ E ' E,M (1) 

to fit the combined ISGRI and JEM-X data, to which we add 
two cyclotron lines (see below). A multiplicative constant is 
applied to allow for the different normalisations of the two in¬ 
struments. Furthermore a systematic error of 2% is applied to 
account for the uncertainties in the response matrices of JEM- 
X and ISGRI. No photoelectric absorption is required at low 
energies, compatible with the fact that we analyse the JEM-X 
data only above 4keV; the addition of a Fe Ka fluorescence 
line does not improve the fits either. 

To fit the cyclotron lines we use Gaussian optical depth pro¬ 
files of the form 

_ l / £ ~ £ cycl ,i 

Gi(E) = Tl ■ e~\ 


(2) 
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Table 2. Fitting parameters. From first to last column: (1) revolution number (see Table0; (2) normalisation factor; (3) powerlaw 
index; (4) folding energy (cutoffpl in XSPEC); (5) energy of the fundamental cyclotron line; (6) width of the line; (7) optical 
depth at the line center; (8) flux in the 2-10 keV band; (9) flux in the 20-40 keV band. 


Rev. 

A 

r 

^fold 

^cycl,l 

cr i 

Tl 

Ecyc\,2 

0-2 

Tl 

2-10 keV 

20-40 keV 




keV 

keV 

keV 


keV 

keV 


erg/cm 2 /s 

erg/cm 2 /s 

273 

n +0.004 

U.4jo _ 0002 

-0 174+ 0005 
' ^-0.004 

7 04 +0 - 01 
' - 0.02 

97 c^+0.04 
z, / . j j_oo4 

4 7o +0 03 
' -0.03 

1.80!°;°! 

51 5 +0 * 3 

9 6 +tu 
^ -0.2 

1 Q 9 + 0.02 

A -^-0.02 

1.89 x 10 ~ 8 

7.31 x 10 ~ 9 

274 

0 35 +0 01 
- 0.00 

-0 26 +0 02 
U,Z,U - 0.02 

6 9 +01 

97 7 + 0.1 
^'-O.l 

4 7 + 0.1 
^ - 0.1 

1-78!^ 

55 0 +L3 
u - 0.2 

11 3 +L3 

11 -0.6 

2 0 +ai 
^ - 0.1 

1.66 x 10- 8 

6.92 x 10 - 9 

278 

0 17+ 001 

U - A ' -0.00 

-0 24 +0 02 

u -^ h -0.07 

7 5 +03 

1 • J -0.3 

28.8!“ 

4 8 +01 

i Q9+0.04 
l ‘ y ^-0.03 

54 6 +L1 
J ^- u -0.8 

q 7+O.7 

^•'-o.s 

2 6 +0 * 3 
z, ‘ - 0.1 

8.53 x 10 ~ 9 

4.37 x 10 “ 9 

284 

0 06 +0 0 ° 
u - uu -o.oo 

-0 4l +000 
U *^ A - 0.01 

7 4 + 0.0 
' - H -o.o 

9 Q a+ 0.0 

Ay ‘ D -o.o 

4 Q+o.o 
^- 0.0 

i 96+0.02 
a -^ u -0.02 

51.0 fix 

7 4 + 0.1 
'•^- 0.1 

2 0 +01 
^• w -o.i 

3.94 x 10 ~ 9 

2.69 x 10 “ 9 

285 

0 04 +0 0 ° 
- 0.00 

-0 53 +0 * 01 

-0.03 

7 9+O.O 
' * -0.3 

9Q 4+0.1 
^ -0.1 

4 4 + 0.0 
^ - 0.1 

1 96+ 004 

1 -0.04 

51.0 fix 

7 q+0.2 

7 -0.2 

2 l +ai 
^ -0.1 

3.15 x 10 ~ 9 

2.41 x 10 ~ 9 

286 

0 03 +0 0 ° 
-0.00 

-0 44 + 005 
U *^ -0.06 

8 -°-a4 

29 5 +0 - 3 

4 5 +01 
^ -0.1 

2 00 +006 

U -0.05 

51.0 fix 

7 Q+0.2 

'■ y -o.s 

9 4 + 0.2 
^ - 0.1 

2.19 x 10 - 9 

1.80 x 10- 9 

287 

0 02 +0 0 ° 
u - uz, -o.oo 

-0 41 +005 
u -^ a -o.or 

83!“ 

29.4!“ 

4 1 + 0.2 

911 +0.08 
11 -0.07 

51.0 fix 

8 . 8 !“ 

9 7+0.3 
' -0.2 

1.51 x 10 - 9 

1.28 x 10- 9 



Fig. 3. JEM-X and ISGRI spectra integrated, from upper to 
lower curves, during revolutions 273,278,284 and 287, respec¬ 
tively. The spectra have been normalised to the Crab spectrum 
of revolution 239 (see text for details). 


where E cyc \ is the line energy, cr the line width, and r the op¬ 
tical depth at the line center (see Kreykenbohm et al. 2Q'(f4 
Pottschmidt et a1. (2005) . The index i identifies the line, i.e i= 1 
for the fundamental at 28 keV and i —2 for the first harmonic at 
50 keV. The second harmonic at 71 keV (Paper I) is not mod¬ 
elled due to insufficient statistics for revolution averaged spec¬ 
tra. The overall spectral model is then given by 

f(E) = F(E ) ■ e~ Gm • e~ GliE) (3) 

where F(E) is given by Eq. Q and G\(E ) and Gi(E)by Eq. 0. 
The same model is used for all observations to keep consistency 
between the parameters from one observation to the next. The 
last observation (revolution 288) is disregarded from this spec¬ 
tral analysis due to poor statistics. 

The main fit parameters are summarised in Table □ The 
2-10keV and 20-40keV fluxes derived from the fitted models 
confirm the two phases set forth in Sect l4. II i.e. first an ex¬ 
ponential decay followed by a linear decrease of the flux. The 
folding energy remains about constant around 7.5 keV, which 
is consistent with the results of Paper I (Table 1) for revolution 
273-274. The power law index f, on the other hand, decreases 



53380 53385 53390 53395 53400 53405 53410 53415 53420 
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Fig. 4. Evolution of the specral fit parameters describing the 
fundamental cyclotron line, a line energy, b width and c depth. 
The horizontal bars on the points indicate the duration of each 
observation. 


from -0.18 in the early revolutions to -0.4 in the later ones, 
in agreement with the hardening observed from the hardness 
ratios (Fig.0. 

The cyclotron line parameters for the fundamental are dis¬ 
played in Fig. El as a function of time, and an example of the 
unfolded spectrum with the fit and residuals for revolution 284 
is shown in Fig. 0 The energy of the fundamental cyclotron 
line increases by 7% from the first to the last INTEGRAL ob¬ 
servations, while the depth of the line increases by 16% and the 
width decreases by 12%. The changes are significant; attempt¬ 
ing to fit the spectrum in revolution 284 with a line at 27.5 keV 
leads to strong residuals around the line (Fig.0). The changes 
in the line parameters cannot be attributed to variations of the 
continuum either; fixing T in revolution 284 to its value of -0.17 
from revolution 273 leads to strong residuals at energies above 
45 keV (Fig-Et)- 

The determination of the parameters describing the 50 keV 
line becomes problematic from revolution 284 on, due to poor 















normalized counts/sec/keV 
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Fig. 5. a Combined spectrum and folded model of data ob¬ 
tained with JEM-X (left) and ISGRI (right) for revolution 284; 
b residuals of the model with r=-0.41 and £f 0 id=7.4 keV, and 
without any Gaussian line applied; c one Gaussian line is in¬ 
cluded at 29.3 keV; d a second Gaussian is added, at 51 keV; e 
same as d, but with the first line fixed at its value of 27.5 keV 
found in revolution 273; f same as d, but with T fixed the value 
of -0.17 found in revolution 273. 

statistics. Therefore, for those revolutions, we fix the line en¬ 
ergy to the value £^ 1,2 =51 keV obtained for Rev. 273 where 
the statistics are the best, and fit the width and depth of the line 
to obtain the best residals. The results are shown in Tabled The 
width and depth of line are seen to be rather stable throughout 
the whole decay of the outburst, around 9keV and 2, respec¬ 
tively. The large error bars on the numbers result from the poor 
statistics. 

4.3. Power spectra 

Several X-ray pulsars and cyclotron line sources are known to 
exhibit low frequency quasi-periodic oscillations (QPOs) (see, 
e.g., Takeshima et al. 1101)1 Takeshima et al. 119941 Heindl et al. 
irm\ . Takeshima et al. ( 1(194 :) first discovered a QPO with a 



Frequency [Hz] 


Fig. 6. 5-20 keV power spectra, averaged over the indicated 
revolutions and normalised using the method described by 
Miyamoto et al. dT99lT) . Poisson noise has not been subtracted 
(see text for more details). The power spectra of revolutions 
278 and 273 have been multiplied by factors of 5 and 25, re¬ 
spectively. Error bars have been omitted for clarity. 


centroid frequency of -0.05 Hz and a relative root mean square 
(rms) amplitude of -5% in 2.3-37.2keV data from the 1989 
outburst of V 0332+53. This feature has recently also been 
reported in the 2-60 keV RXTE -PCA power spectral density 
functions (PSDs) obtained during the 2004 outburst decay (Qu 
et al. 120051 . Furthermore Qu et al. found a QPO at -0.22 Hz, 
i.e. centered on the fundamental pulsar frequency, in addition 
to the narrow peak at that frequency. The rms strength of this 
QPO also amounts to a few per cent. In the following we show 
that the JEM-X and ISGRI observations confirm the detection 
of both features during the outburst decay, and we analyse their 
changes with time. Note that, to our knowledge, these results 
represent the first QPOs reported from INTEGRAL JEM-X 
and ISGRI data in the literature. 

The background subtracted 5-20 keV and 20—40keV 
lightcurves of each revolution were rebinned to a time resolu¬ 
tion of 0.5 s and power spectra were produced from 512 s long 
lightcurve segments using Fast Fourier Transforms. Averaged 
power spectra for selected revolutions are shown in Figs. [6] 
and [7] 

For coded mask instruments the determination of the white 
noise level as well as the influence of the background flux on 
the power spectrum can be complex. The white noise compo¬ 
nent in JEM-X lightcurves, e.g., is known to be non-Poissonian. 
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Fig. 7. 20-40 keV power spectra, averaged over the indicated 
revolutions. The spectra are arbitrarily shifted in the y-axis for 
easy comparison. Note that the ISGRI power spectra of revolu¬ 
tion 278 are not shown due to solar flaring giving rise to many 
gaps in the lightcurve after GTI filtering. 

This is also true for the ISGRI which consists of 16000 pix¬ 
els, each having different noise charactersistics. Using the nor¬ 
malization described by Leahy et al. l fM5l we find the noise 
level to be about one order of magnitude higher than expected 
from Poisson noise. However, source features are well appar¬ 
ent above the noise level, e.g., the 0.05 Hz QPO is clearly vis¬ 
ible in revolution 284. Therefore, and since the broad-band 
PSD (which is not to be discussed here) is dominated by the 
noise component up to at least 0.1 Hz, we did not attempt to 
correct the power spectra for white noise. The JEM-X power 
spectra shown in Fig.[6]are normalised according to Miyamoto 
et al. ( llQOll and a comparison of the strength of features 
above the noise level between the JEM-X PSDs and the quasi- 
simultaneous PCA PSDs of Qu et al. 12005) finds comparable 
values 1 . 

Strong changes in the shape of the power spectrum can be 
seen over the outburst decay, where the varying amplitude ra¬ 
tios at the pulsar frequency and its harmonics reflect a strong 
change of the pulse profile (Zhang et al. B005i l. The PSDs 
shown in Figs. [6] and 0 represent the main stages of the evo¬ 
lution. Although, for the reasons discussed earlier, the white 
noise level is not accurately estimated in the INTEGRAL data, 
the JEM-X and ISGRI PSDs for revolutions 0273 and 0284 do 
compare well and show similar features as a function of fre¬ 
quency. 

In revolution 273 (Jan. 8-9) the power spectrum is domi¬ 
nated by signatures of the pulsation frequency, especially a very 
strong first harmonic. In addition weak red noise is present. The 
peak at the fundamental pulsation frequency is slightly broad¬ 
ened at the bottom, especially in the ISGRI PSD. However, 
the 0.22 Hz QPO is not as obvious as a few days later in the 

1 This can only be a very rough comparison since different energy 
bands are involved and since the power spectrum of V 0332+53 obvi¬ 
ously changes rapidly during the decay. 


Jan. 15-19 PCA power spectrum of Qu et al. (120051 . In rev¬ 
olution 278 (Jan. 23-24) the fundamental pulsation at 0.22 Hz 
is not seen anymore but a weak broad QPO feature is present 
at this frequency. The first harmonic is still clearly visible. The 
strength of the red noise component more than doubled at low 
frequencies and there is an indication of additional power in 
the frequency region of the 0.05 Hz QPO. Principally the over¬ 
all PSD shape is very similar to the Jan. 15-19 PCA power 
spectrum mentioned above. Finally, in revolution 284 (Feb. 9- 
10) the 0.05 Hz QPO is clearly visible, especially in the JEM-X 
power spectrum. The red noise is still strong while the pulsa¬ 
tion fundamental is back as a narrow feature, stronger than the 
first harmonic. Fitting the low frequency QPO in revolution 284 
with a Lorentzian gives (0.049 + 0.001)Hz, with a FWHM of 
(0.023 + 0.005) Hz for JEM-X and (0.049 + 0.003) Hz with a 
FWHM of (0.021 + 0.008) Hz for ISGRI, where the underlying 
continuum can be modelled with a flat power law component 
and a very low frequency Lorentzian. 

5. Discussion 

V 0332+53 is a highly magnetised pulsar which has a magnetic 
field of about 2.7 x 10 12 G (Paper I, Pottschmidt et al. 2005 1 
and orbits a massive 08-9Ve companion with a period of about 
34 days. The X-ray outbursts from this system are observed 
with a time recurrence of about 10 years. They are believed 
to be triggered by a sudden increase of the mass ejection rate 
from the massive star, the signature of which has been observed 
in the optical brightening of the system in 2004 (Goranskij & 
Barsukova l2004l i. The ejected mass feeds an accretion disk 
around the pulsar, from which it is funneled by the magnetic 
field towards the polar caps and falls on the surface of the neu¬ 
tron star. 

The analysis of the 2005 outburst presented in Sect. □p ro- 
vides new hints towards a better understanding of V 0332+53 
within this basic picture. The X-ray lightcurve (Sect. l4~fl dis¬ 
plays a characteristic exponential decay with a time scale of 20 
to 30 days, followed by a linear decrease. This feature may pro¬ 
vide some information on the disk properties. The X-ray spec¬ 
trum (Sect. IP is well fitted by a folding power law modified 
by two cyclotron lines. The 27 keV energy of the fundamental 
cyclotron line is shown to increase by 7% during the outburst 
decay and to become narrower by 12%. At the same time, its 
depth increases by 16%. This information would provide some 
insight into the resonance scattering region near the polar caps. 
Finally, the appearance of the 0.05 Hz QPO only from the mid¬ 
dle of the decay on (Sect. IP provides further information on 
the properties of the matter around the neutron star. Each of 
those observational facts are successively discussed in the fol¬ 
lowing sections. 

5.1. Lightcurve 

An exponential decay of the flux is observed in dwarf novae 
and in soft X-ray transients on timescales from a few days for 
the former to several weeks for the latter (e.g. Lasota 119961 
Tanaka & Shibazaki 1 996 1 King & Ritter ( 1998 ) show how 
illuminated disks in low-mass X-ray binary stars can produce 
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such an exponential decay on the time scales observed for the 
soft X-ray transients. In that picture, the central star illuminates 
the disk and keeps its temperature high. The disk empties pro¬ 
portionally to its mass, leading to an exponential decay of the 
emitted flux. They further show that the lightcurve would be 
characterised by a linear decrease at later times. The change 
to a linear regime is triggered by a change in the conditions in 
the disk. In low-mass X-ray binaries the change corresponds to 
the decrease of the disk temperature below the temperature of 
partial H ionisation, such that the high temperature cannot hold 
over the entire disk anymore. This scenario of X-ray irradiated 
accretion disks has subsequently been confirmed by numerical 
calculations (Dubus et al. n999ii2orm> . 

It is tempting to make the parallel with the results found for 
V 0332+53. The source displays an exponentially decreasing 
flux followed by a linear decrease, similar to the disk-irradiated 
systems. This is the first time, to our knowledge, that such a 
behaviour is observed in a high-mass X-ray binary. Since the 
X-ray flux is proportional to the mass accretion rate, this picture 
suggests that the accretion rate is proportional to the mass of 
the disk during the exponential decay phase. The transition to 
a linear decrease phase would be triggered by an as yet un¬ 
identified change of conditions in the disk. 

The two different exponential decay times observed in the 
lightcurves of V 0332+53, of 30 days above 20keV and 20 
days below (Sect. EED suggest the presence of at least two re¬ 
gions contributing to the X-ray continuum of that source. The 
harder component would originate from the accretion column 
near the resonant scattering region, while the softer component 
would originate from a region located, for example, higher in 
the accretion column. The observed increase of the depth of 
the cyclotron line with time is also compatible with this sce¬ 
nario; the softer component decays more rapidly and hence its 
tail would contribute relatively less to the high energy part of 
the spectrum from which the photons are scattered by the elec¬ 
trons in the magnetic field. Basically, one time scale governs 
the rate of mass flow onto the neutron star through the disk, 
leading to the exponential decay followed by a linear decrease 
(see above); another time scale governs the spectral changes of 
the emission that may come from different heights within the 
accretion column. 

5.2. Spectrum 

Time evolution of cyclotron resonance line parameters has al¬ 
ready been reported in the literature for some binary systems. 
Mihara et al. ( 1998 ;) report a relation between the line energy 
and the X-ray luminosity for 5 pulsars observed by Ginga. The 
resonance energy of 4U 0115+53, for example, increased by 
40% as the luminosity decreased by a factor of 6 between 1990 
and 1991. The 1990 outburst of V 0332+53 is also reported 
in that paper, with a 10% decrease of the line energy from 
30.0+0.5 keV to 27.2+0.3 keV as the Ginga count rate in the 
3-37 keV band increased by a factor of 1.6. These authors at¬ 
tribute the change of the line energy to that of the height of 
the scattering region in a dipole magnetic field, leading to a 
height change of ~l.lkm for 4U 0115+53 and ~330m for 


V 0332+53. Kreykenbohm et al. (12004I) . on the other hand, re¬ 
port a variation of the cyclotron energy of GX 301-2 with the 
pulse phase of the pulsar, and attribute the change to different 
viewing angles of the accretion column where the line origi¬ 
nates. 

In the case of the 2005 outburst decay of V 0332+53 anal¬ 
ysed in this paper, the increase with time of the energy and 
depth of the cyclotron fundamental line, together with the de¬ 
crease of its width can also be understood as resulting from 
the variation of the scattering region characteristics in a dipole 
magnetic field above the polar caps of the pulsar. It is inter¬ 
esting to note that Pottschmidt et al. d2005b . analyzing the 
RXTE data of V 0332+53 taken two weeks before the first 
INTEGRAL observation, derive Ei = 25.2 keV for 2004 
December 24-26, which confirms the trend of increasing line 
energy with time found in this paper. Let us thus assume that 
the region where the cyclotron resonance scattering takes place 
extends over some fraction of the accretion column. Then the 
region would cover some range of magnetic strengths, and a 
broad line shape would result from the superposition of a suc¬ 
cession of narrow lines, each originating from a different height 
and magnetic field strength in the accretion column. As the ac¬ 
cretion rate and hence the flux decreases with time, the exten¬ 
sion of the resonance scattering region is expected to decrease 
along with the total width of the line. At the end of the de¬ 
cline only the region close to the neutron star surface, where 
the magnetic field is stronger, would contribute to the reso¬ 
nance scattering; this would lead to a narrower line shape at 
a higher energy. The observed energy Ei indeed increases by 
17% from the RXTE observation of December 2004 to the last 
INTEGRAL observation in February 2005. If we consider a 
dipole approximation for the magnetic field, we have B ~ r ! 
in a non-relativistic description, where r is the radius to the 
center of the neutron star. The extent of the resonance scatter¬ 
ing region at the peak of the outburst would then, in this first 
scenario, be about 6% of the stellar radius, i.e. about 500 me¬ 
ters for typical radii of neutron stars. A change in the scattering 
energy between two different observations of Her X-l, from 
34keV to 41 keV, has been discussed by Gruber et al. d2?Xm> : 
they conclude that the range over which the scattering occurs 
is of the order of 1 km in height, if their observed line width of 
5 keV is largely due to a mix of magnetic field strengths in a 
dipole field. In the case of V 0332+53, we additionally observe 
an evolution of the width and depth of the line with time. 

Alternatively, we can argue that the fundamental line orig¬ 
inates from a unique narrow region in the column density, but 
that the line shape is quite complex and variable with time, re¬ 
sulting in variations of the line width and depth. An increase 
in the line energy would then be related to a displacement with 
time of the scattering region. In the case of V 0332+53, the 
observed 17% increase of E! would imply a move of the lo¬ 
cation of the scattering region by 500 meters in the accretion 
column down to the surface of the neutron star. However, in 
this second scenario, the line profile is expected to remain ap¬ 
proximately constant with time, which is not supported by the 
INTEGRAL observations. Also, if the scattering region gets 
closer to the neutron star, higher order terms in the B-field may 
increase the line width, which is opposite to the trend observed 










8 


N. Mowlavi, I. Kreykenbohm, S. E. Shaw et al.: V 0332+53 during the 2005 outburst decline 


by INTEGRAL. We therefore would rather support the first 
scenario. 

5.3. Quasi periodic oscillations 

The 0.22 Hz QPO observed around the spin frequency in both 
the INTEGRAL (Sect. [H} and RXTE (Qu et al. l2(T05t data 
must originate from a region around the neutron star co-rotating 
with the pulsar. This region can either be located at the inner 
radius of the accretion disk near the Alfven radius, or in the 
corona above the hot spots, or somewhere in between along the 
magnetic field through which the matter is tunneled. 

The origin of the lower frequency QPO at 0.05 Hz, on the 
other hand, is more debatable. Low frequency QPO-type vari¬ 
abilities are known in several X-ray pulsars (see van der Klis 
120041 for a review). Although there is no general consensus on 
the origin of QPOs, most models do invoke the presence of 
some sort of inhomogeneities in the accretion disk. With a typ¬ 
ical neutron star mass Mns = 1.4M 0 and a spin of Pns —4s 
for V 0332+53, the co-rotation radius where the angular veloc¬ 
ity of the magnetosphere and the Keplerian velocity of the disk 
are equal is 

flcorot- 1500 (—) ' km - 4230km. (4) 

On the other hand, the radius at which the Keplerian frequency 
is equal to that of the 0.05 Hz QPO is ~12350km. The scenario 
of having the 0.05 Hz QPO somehow related to the part of the 
disk at this distance, i.e. three times further away than the co¬ 
rotation radius, is thus plausible. The temperature of the disk 
cannot be responsible for the observed emission, since no black 
body emission is found in the spectrum. The QPO could result 
from an occultation of the beam by intervening matter from the 
disk at that distance, as suggested by Heindl et al. (US to 
explain a 2 mHz QPO observed in 4U 0115+63. 
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